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Abstract 


Scattered  light  from  particles  near  a  plane  interface  has  important  impli¬ 
cations  for  many  branches  of  science.  Applications  range  from  radar  de¬ 
tection  and  remote  sensing  to  clean-room  monitoring  and  quality  control 
in  the  manufacture  of  silicon  wafers.  The  most  fundamental  system  of  this 
type  is  that  of  a  sphere  near  a  plane  interface.  Calculating  the  scattering 
is  complicated  because  the  boundary  conditions  at  the  sphere  and  plane 
surfaces  must  be  satisfied  simultaneously,  and  these  two  systems  represent 
two  fundamentally  different  geometries.  In  this  report  we  present  an  an¬ 
alytical  solution  to  this  problem  that  retains  the  physics  and  may  be  used 
as  a  basis  for  more  numerically  intensive  techniques  necessary  when  the 
system  constituents  become  irregular. 
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1.  Scattering  System 


Physically,  we  can  consider  four  field  components  that  interact  with  the 
system  of  a  sphere  near  a  plane  interface  shown  in  figure  1.  For  one  compo¬ 
nent,  the  incident  light  strikes  the  sphere  and  is  scattered  into  the  far-field 
region  beneath  the  plane  interface.  For  a  second  component,  the  incident 
light  reflects  off  the  plane  interface  before  it  is  scattered  by  the  sphere  into 
the  far-field  region  beneath  the  plane  interface.  For  a  third  component,  the 
incident  light  strikes  the  sphere  and  is  scattered  toward  the  plane  interface, 
which  reflects  this  light  toward  the  far-field  region  beneath  the  plane  inter¬ 
face.  And  for  a  fourth  component,  the  incident  light  reflects  off  the  plane 
interface  before  it  is  scattered  by  the  sphere  toward  the  plane  interface;  the 
plane  interface  reflects  this  light  toward  the  far-field  region  beneath  the 
plane  interface.  The  total  scattered  field  can  be  found  as  the  superposition 
of  these  four  components,  since  it  is  widely  known  how  a  plane  wave  in¬ 
teracts  with  isolated  spheres  and  plane  interfaces.  However,  this  "double 
interaction"  model  does  not  include  the  interaction  fields,  i.e.,  fields  that 
are  scattered  by  the  sphere,  reflected  off  the  plane  interface,  and  are  once 
again  incident  upon  the  sphere.  To  include  the  effects  of  this  interaction,  we 
must  use  more  complicated  models. 

In  this  report,  we  discuss  a  multipole  solution,  based  on  image  theory,  to 
the  problem  of  a  sphere  in  proximity  to  a  plane  interface.  This  type  of  solu¬ 
tion  is  valuable,  since  it  retains  the  individual  field  components  interacting 
with  the  individual  system  components  that  many  numerical  techniques 
cannot.  In  this  way,  the  physics  is  preserved.  The  incident,  scattered,  and  in¬ 
teraction  fields  are  expanded  in  terms  of  vector  spherical  harmonics  about 
coordinate  systems  centered  on  the  sphere  and  its  image.  The  boundary 
conditions  on  each  subsystem  can  then  be  satisfied  individually.  Because 
the  field  expansions  are  in  terms  of  vector  spherical  harmonics,  it  is  rela¬ 
tively  easy  to  satisfy  the  boundary  conditions  on  the  spherical  interface.  To 
satisfy  the  boundary  conditions  at  the  plane  interface,  an  image  (or  interac¬ 
tion)  field  is  introduced  that  is  physically  equivalent  to  the  scattered  field 
from  the  sphere  that  reflects  off  the  plane  interface  and  is  incident  upon  the 
sphere.  For  a  perfectly  conducting  plane  interface,  this  interaction  field  is 
the  image  of  the  scattered  field  from  the  sphere.  When  the  half-space  not 
encompassing  the  particle  is  dielectric  or  lossy,  some  approximation  or  nu¬ 
merical  technique  must  be  used  to  represent  this  interaction  field. 
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Figure  1.  Diagram  of 
scattering  system.  A 
sphere  of  radius  rgph , 
complex  refractive  index 
rriaph,  lies  a  distance  d 
beneath  a  plane  interface 
separating  a  medium  of 
refractive  index  mi  below 
interface  {zi  <  d)  from  a 
region  of  complex 
refractive  index  m2  above 
interface  {zi  >  d).  A  plane 
wave,  whose  wave  vector 
lies  in  the  xi  —  21  plane 
and  is  oriented  at  angle  9o 
with  respect  to  21  axis,  is 
incident  upon  system.  To 
help  satisfy  boundary 
conditions,  an  image 
coordinate  system 
{x2,y2,  Z2)  is  introduced  a 
distance  2d  above 
coordinate  system 
centered  on  sphere. 
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2.  Field  Expansions 


Figure  1  shows  the  scattering  system.  The  isotropic,  homogeneous  sphere 
of  refractive  index  msph  and  radius  rsph  is  located  in  air,  centered  on  the 
(xi,  yi,  zi)  coordinate  system  located  a  distance  d beneath  a  plane  interface 
separating  region  2  of  complex  refractive  index  m2  {zi  >  d)  from  region 
1  of  real  refractive  index  mi  {zi  <  d).  We  denote  ki,  k2,  and  fcsph  to  be  the 
propagation  constants  of  waves  in  region  1,  region  2,  and  within  the  sphere. 
To  satisfy  the  boundary  conditions  at  the  plane  interface,  we  introduce  an 
image  coordinate  system  (x2,  y2,Z2),  located  a  distance  d  above  the  plane 
interface  on  the  zi  axis.  This  coordinate  system  is  the  source  of  the  image 
or  interaction  fields.  Incident  upon  the  scattering  system  is  a  plane  wave  of 
wavelength  A  whose  wavevector  ko  is  oriented  in  the  xi  —  zi  plane  at  angle 
9o  with  respect  to  the  zi  axis.  The  electromagnetic  fields  may  be  expanded 
in  terms  of  the  vector  spherical  harmonics  centered  about  the  j  coordinate 
system  {j  =  1,2): 


M^^l{krj)  =  Oj  z!^P\krj)iTTmn{cos9j)exp{im<j)j) 

-  zlf'> {krj)fmn  (cos  9j)  exp  (imipj) 


=  Tj  -^z!^P\krj)n{n  +  l)P^{cos9j)exp{im(l)j) 


+  %  I  rjZ^\krj)  Tran  (cos  9j)  exp  (im(/)j)  | 

+  I  ^  ^  {krj)  TTmn  (cos  9j)  exp  |  , 

where  Zn\krj)  are  the  spherical  Bessel  functions  of  the  first,  second,  third, 
or  fourth  kind  {p  =  1,2, 3, 4),  and 

777 

TTmn  (cOS  9)  =  ^—P^{cOS  9)  (3) 

Sin  6 


fmn  (cos  9)  =  —  P™(C0S  9) 
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(5) 


P-(cos0)  = 


'(2n  +  l)(n  —  m)\ 


2{n  +  m) 


pm 

■^n 


(cos  9), 


where  P^{cos6)  are  the  associated  Legendre  polynomials.  We  assume  a 
time  dependence  of  exp(— iwt)  throughout. 


To  satisfy  the  boundary  conditions,  we  use  the  linearity  property  of  the 
electromagnetic  fields  and  consider  the  components  of  the  electromagnetic 
field  individually.  These  components  are  each  expressed  as  a  multipole  ex¬ 
pansion  of  vector  spherical  harmonics,  and  the  boundary  conditions  at  the 
interfaces  of  the  sphere  and  of  the  plane  are  satisfied  simultaneously.  We 
first  consider  the  incident  fields. 
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3.  Incident  Field 


The  incident  electric  field  that  strikes  the  sphere  is  expanded  in  vector 
spherical  harmonics  about  the  sphere  coordinate  system  as 

cxD  n 

Einc  =  E  E  (fcirO+a^NW  (6) 

r^  =  ^  m  .=—r). 


where  a^n  are  the  incident  field  coefficients.  This  component  of  the  electric 
field  is  made  up  of  two  parts:  (1)  the  plane  wave  that  strikes  the  sphere 
directly  and  (2)  the  plane  wave  that  strikes  the  sphere  after  it  is  reflected 
off  the  plane  interface.  We  consider  these  subcomponents  individually,  and 
because  we  are  dealing  with  a  linear  system,  the  incident  field  coefficients 
are  each  composed  of  two  elements: 


„(i)  ^  „(i)o 

^ran  ^ran 


where  a^n  are  the  coefficients  of  the  incident  plane  wave  that  strike  the 
sphere  directly,  and  am^  are  the  coefficients  of  the  incident  plane  wave  that 
were  reflected  off  the  plane  substrate.  In  general,  we  consider  two  different 
polarization  states  of  the  incident  field:  (1)  A  unit-normalized  TE-incident 
(transverse  electric)  plane  wave  of  the  form  =  expi/ci(zi  cos  0o  + 
xisin^o^yi  and  (2)  a  unit-normalized  TM-incident  (transverse  magnetic) 
plane  wave  of  the  form  =  exp \ki{zi  cos  9o  +  xi  sin  do)  x  (xi  cos  9o  — 
zi  sin  9o)-  We  denote  the  coefficients  for  these  two  cases  with  an  additional 
superscript:  Omn and  amn’^^ ■  We  can  find  these  coefficients  using  the 
orthogonality  properties  of  the  vector  spherical  wave  functions.  The  first 
coefficient  can  be  expressed  as 


ri)o,TE  ^  /ETE.MW*(fcir-i)dE 

/M(^^(feiri)  •MW*(A:iri)dE’ 

where  the  asterisk  denotes  the  complex  conjugate  and  the  integral  is  taken 
over  all  space.  With  a  little  bit  of  algebra  and  calculus,  we  find 


,(l)oTE  _ 


n{n  +  l) 


,(cos  9o). 
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-2i'^ 


Similarly,  the  coefficients 


,(2)o,TE  _ 


/ETE.N^];(fciri)dy  ^  -2i" 

n{n  +  l) 


TTr 


I  (cos  9o 


(10) 


It  is  a  simple  exercise  of  Maxwell's  equations  to  show  that 


and 


„(l)o,TM  _  :  (2)o,TE 
^mn  ^mn 

(11) 

„(2)o,TM  _  ■  (l)o,TE 
^mn 

(12) 

Finally,  we  must  determine  the  reflected  components  of  the  incident  field 
coefficients.  The  solution  to  the  scatter  of  a  plane  wave  by  an  infinite  plane 
interface  is  a  classic  problem  of  electromagnetism  solved  by  Fresnel  in  the 
nineteenth  century  The  amplitude  and  phase  of  the  reflected  plane  wave 
are  dependent  on  its  polarization  state,  incident  angle,  and  the  refractive 
index  of  the  media  at  either  side  of  the  interface.  The  Fresnel  reflection  co¬ 
efficients  for  a  plane  wave  incident  at  angle  a  on  the  interface,  traveling 
from  medium  i  and  reflecting  off  the  interface  at  medium  j,  are 


(13) 


where  the  superscripts  on  the  reflection  coefficients  denote  the  polarization 
state  of  the  incident  field.  It  should  be  noted  that  these  particular  sets  of 
equations  are  valid  when  the  reflected  electric  field  component  is  pointed 
in  the  same  direction  as  the  incident  electric  field  component  for  both  TE 


and  TM  polarization 

states  when  a  =  0.  From  the  relation 

(15) 

we  find 

TTmn  {-x)  =  (-!)'"+”'  TTmn  (x) 

(16) 

and 

fmn  {-x)  =  fmn  {x)  ■ 

(17) 

Using  these  expressions,  we  can  write  the  coefficients  of  the  reflected  field 
in  terms  of  the  coefficients  of  the  direct  field: 


(-1)^+^+”^  Rjf  {Oo)  exp  (2i/cidcos  6»o)  (18) 
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for  the  TE-incident  plane  wave,  and 


„(i)TM  ^  (i)o,TM 
^mn  ^mn 


(-1) 


j+n+m  ^TM  {2\kidcOS  9o 


(19) 


for  the  TM-incident  plane  wave.  The  total  incident  field  coefficients  can  be 
written  as 


„0)TE  ^  (i)o,TE 
^mn  ^mn 


1  +  (—1)-?+"^+™'  (0q)  exp  (2i/ci(icos  6o) 


(20) 


for  the  TE-incident  plane  wave,  and 


a 


(i)TM 

mn 


=  a 


(i)o,TM 

mn 


1  -|-  (— i?™  (6*o)  exp  (2i/ci(icos  9o) 


(21) 


for  the  TM-incident  plane  wave.  Equations  (20)  and  (21)  are  the  vector 
spherical  harmonic  expansion  coefficients  for  the  total  field  of  a  plane  wave 
reflecting  off  a  plane  interface,  and  are  nothing  more  than  Eresnel  theory  ex¬ 
pressed  in  a  spherical  coordinate  system.  The  phase  term  exp  (2i/cidcos  9o) 
exists  in  these  expressions  because  the  coordinate  system  from  which  the 
waves  are  expanded  is  located  a  distance  d  beneath  the  plane  interface. 
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4.  Fields  at  the  Spherical  Interface 


The  next  step  in  the  derivation  is  to  satisfy  the  boundary  conditions  at  the 
spherical  interface.  All  the  components  of  the  electromagnetic  field  must  be 
expanded  in  vector  spherical  harmonics  as  we  did  in  the  previous  section 
for  the  incident  field.  The  scattered  field  component  is  expanded  as 

oo  n 

Esca=E  E  (fciri)  +  6(^),NW  (fciri).  (22) 

72=1  m=—n 

(o\ 

For  the  chosen  time  dependence,  the  spherical  Hankel  functions  zn  {km)  = 
hn  ^  (kivi)  represent  outgoing  spherical  waves  and  satisfy  the  boundary 
conditions  at  infinity  Because  there  is  a  pole  at  the  origin,  this  expansion 
can  only  be  used  in  the  region  exterior  to  the  spherical  interface. 

The  electric  fields  inside  the  sphere  volume  (ri  <  rsph)  are  expanded  in 
vector  spherical  harmonics  about  the  sphere  coordinate  system  as 

OO  n 

Esph  =  E  E  c^MU^sphri)+c(^iN«  (fcsphri),  (23) 

n=l m=—n 


where  Cmn  are  the  internal  field  coefficients.  The  radial  dependence  of  the 
vector  spherical  harmonic  expansion  must  be  in  terms  of  Bessel  functions 

Zn  ^  {ksphn)  =  jn  {ksphn)  . 

Finally,  an  image  or  interaction  field  must  be  introduced  to  help  satisfy  the 
boundary  conditions  at  the  plane  interface.  This  field  is  expanded  in  a  series 
of  outgoing  spherical  waves  centered  on  the  (x2,  y2,  ^2)  coordinate  system: 

OO  n 

Eima  =  E  E  (feir2)  +  d(^MUfcir2).  (24) 

n=l  m=—n 

Physically,  this  field  corresponds  to  the  interaction  of  the  sphere  with  the 
plane  interface:  the  scattered  field  from  the  sphere  reflects  off  the  plane  in¬ 
terface  and  interacts  with  the  sphere  again.  To  satisfy  the  boundary  condi¬ 
tions  on  the  spherical  interface,  it  is  convenient  to  expand  this  image  field 
on  the  (xi,  2/1,  zi)  coordinate  system.  Such  a  transformation  of  coordinate 
systems  can  be  performed,  since  the  vector  spherical  harmonics  represent 
a  complete  orthogonal  expansion.  We  discuss  details  of  this  expansion  and 
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satisfying  the  boundary  conditions  on  the  plane  interface  in  the  next  sec¬ 
tions.  Because  there  is  no  pole  in  the  interaction  field  when  ri  =  0,  the 
radial  dependence  of  the  vector  spherical  harmonic  expansion  must  be  in 
terms  of  Bessel  functions  Zn'^  {kiri)  =  (fciri): 

oo  n 

Ei„,a=E  E  (fciri)  +  eEN«  (25) 

n=l m=—n 

Satisfying  the  boundary  conditions  at  the  spherical  interface  is  a  relatively 
straightforward  process.  This  has  been  done  well  in  so  many  texts  that  it 
would  be  tedious  to  repeat  the  process.  Since  many  applications  call  for 
scattering  from  nonspherical  particles,  it  is  worthwhile  to  outline  a  more 
general  formulation  of  this  scattering  solution.  The  solution  when  the  par¬ 
ticle  is  spherical  will  then  be  provided.  We  express  the  fields  incident  upon 
and  scattered  by  the  particle  in  vector  spherical  harmonics.  Any  field  inci¬ 
dent  upon  a  three-dimensional  particle 

Einc  =  pW  mW  (A;,ri)  +  (fein)  (26) 


has  a  corresponding  scattered  field,  which  can  be  expanded  in  terms  of 
vector  spherical  harmonics 


OO  n 


Esca  =  E  E  +Pmn2)^nm'n']  (^iri)  (27) 

n'=l  m'=—n' 

-L  L(i)  7"21  (2)  j.22  ,1  iA-iriV 

'  ranm'n'  '  Jrran-^mnra'n' \  m'n'  ’ 


that  is,  the  scattered  field  coefficients  bmn  are  related  to  the  incident  field 
coefficients  via  some  transition  matrix  (or  T-matrix)  ’ 


OO  n' 

=  E  E 


„(1)  „(2)  j^j2 

irmn-^mnm'n'  '  mnm'n' ' 


n'=l  m'=—n' 


(28) 


In  a  multiple  system  like  that  of  a  sphere  and  a  plane  interface,  the  total 
field  incident  upon  the  particle  is  composed  of  the  incident  plane  wave 
and  the  interaction  field,  and  the  scattered  field  coefficients  are  expressed 
as 

OO  n'  OO  n 

f-a  =  E  E  EE 

n'=l  m'=—n'  n=l  m-—n 


^mn  '  ^mn 


mnm'n'  ' 


a(2)  +  e(2) 
^mn  '  ^mn 


nj‘^ 


(29) 


This  equation  expresses  the  scattering  coefficients  in  terms  of  the  T-matrix 
of  the  isolated  particle.  This  is  convenient,  since  T-matrix  routines  are  read¬ 
ily  available  for  many  particle  systems.  For  certain  types  of  particles,  the 
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T-matrix  can  be  simplified.  The  T-matrix  of  axially  symmetric  particles 
obeys  radially  symmetric  particles,  = 

Snn'Smm'T^nm'n'^  and  there  is  no  mode  mixing,  i.e.,  =  T^nm'n'  =  0- 

Spheres  are  a  subgroup  of  this  latter  particle  type,  and  the  T-matrix  is  sim¬ 
ply  composed  of  the  Mie  scattering  coefficients  /h 


1301 

ksph'4’n  (fephrsph)  (^l^sph)  (^l^sph)  V’n  (fephrsph) 

^sphV’n  (^sph^sph)  (^l^sph)  (^l^sph)  (^sph^sph) 


fn  (31) 

^iV’n  (^sph^sph)  V’n  (^l^sph)  fephV^n  (^l^sph)  V’n  (fephrsph) 
(^sph^sph)  ?n  (^l^sph)  fcsphV^n  (^l^sph)  (^sph^sph) 

where  ipn  (x)  =  xjn  {x)  and  {x)  =  xhn'^  (x)  are  the  Riccati-Bessel  func¬ 
tions  and  the  primes  denote  derivatives  with  respect  to  the  argument.  The 
other  terms  of  the  T-matrix  are  zero. 


and 

^  mnmn 


10 


5.  Translation  Addition  Theorem 


The  translation  addition  theorem  is  used  to  translate  a  set  of  vector  spheri¬ 
cal  harmonics  from  one  coordinate  system  onto  another  coordinate  system. 
We  need  to  find  in  this  derivation  a  relationship  between  the  two  sets  of  in¬ 
teraction  coefficients  given  by  equations  (24)  and  (25).  The  fields  expressed 
by  these  equations  are  identical,  but  they  are  expanded  about  different  co¬ 
ordinate  systems.  For  a  translation  along  the  negative  zi  axis  a  distance  2d 
with  no  rotation,  the  vector  spherical  harmonics  are  related  by 

OO 

m2  (km)  =  E  a2’”>m2  (fcin)  +  (tin)  (32) 


n2  (km)  =  E  Bi”’”’M2-  (fcii-l)  +  2”’"’n2-  (fcii-l) .  (33) 

n'=l 

This  relationship  is  valid  in  the  region  where  ri  <  \2d\  ,  i.e.,  out  to  the  pole 
located  at  r2  =  0.  We  can  use  these  equations  to  find  the  relationship  be¬ 
tween  the  dmn  and  Cmn  coefficients.  Substituting  equations  (32)  and  (33) 
into  equation  (24)  and  comparing  with  equation  (25)  yields 


(_xj 

(1)  _  j(l)  Aim,n')  j(2)  o{m,n') 


p(2)  _  j(l)  I  ,(2)  A{m,n') 

^mn  /  V  *^mn'  n  '  ^mn'  n  5 

n'=l 

and  for  a  spherical  particle,  equation  (29)  can  now  be  expressed  as 


hU)  =  V  V  V  4(™X)  I  j(2)  n(mX)l  1  .1 

^mn  /  .  /  .  J  “mn  '  /  .  '  ^mn'^n  )[■'”' 

n=l  m=—n  K  n'=0  ) 

+  j o2  +  £  KwBi’”-"''  +  iTA'.™'”'’]  I  fn  ■  (36) 

I  n'=0  ) 

It  is  beyond  the  scope  of  this  report  to  provide  a  complete  derivation  of 
the  vector  translation  coefficients  and  ;  however,  it  is  certainly 
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worthwhile  to  write  out  a  set  of  relations  from  which  these  may  be  calcu¬ 
lated.  The  vector  translation  coefficients  can  be  expressed  in  terms  of  the 
scalar  translation  coefficients  that  are  used  to  translate  the  spherical 

harmonics  from  one  coordinate  system  onto  another  coordinate  system: 

CXD 

h!^\kr2)Pn  {cos  62)  exp  {im(p2)  =  (cos  ^i)  exp  {im(j)i)  .  (37) 

n'=0 


The  scalar  translation  coefficients  can  be  found  via  various  recursion  rela¬ 
tions.  The  starting  point  for  calculating  the  scalar  coefficients  is 


which  results  from  the  expression 

00 

^o^^(A:ir2)  =  {‘^n  +1)  jn'{kiri)P^\2kid). 

From  here,  the  relation 


n'=0 


^(0,n+l)  _ 


1 


2n  -h  3  I  ,  2n  +  l  „(o,n) 


{n  +  l)\l  2n'  +  1 
-|-  n 


n 


2n'  -  1 


(38) 


(39) 


(40) 


l2n'  1  (o,n-i)  _  ,  ,  y  j  2n  +  1~  (o,n) 

V2n-1  ^  +  ^V2n'  +  3  -'+1 


can  be  used  to  calculate  the  elements  .  The  remaining  terms  in  the 
matrix  can  be  calculated  using  the  relation 


^(m,n)  _  jin'  -m  +  l){n'  m){2n'  1)  (^rn-i,n) 

y  (n  —  m -|- l)(n -|- m)(2n' -|- 1) 

_2k  rl  !  {n' -m  +  2){n' -m  +  1)  ~ 

^  \l  {n  —  m  +  l)(n  -|-  m){2n'  -|-  l)(2n'  -|-  3) 


(41) 


—  2kid^ 

It  should  be  noted  that 


(n'  -|-  m)  (n'  -|-  m  —  1) 


{n  —  m  +  l)(n  -|-  m){2n'  -|-  l)(2n'  —  1) 


C'n'-l 


yj{m,n)  _  yj(-m,n) 


(42) 


We  can  find  the  vector  translation  coefficients  using 


2kid  /(n'  —  m  +  l)(n'  +  m  +  1)  ^(m,n) 


_  Qyin.iL)  _ 


n‘ 


_  _ 

+  l\l  {2n' +  l){2n' +  3) 


2kid  I  {n' -  m){n' +  m)  ^(rn,n) 


n 


(2n'  -h  l)(2n'  -  1) 


(43) 
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j-){m,n)  _  ‘2\k\Tnd 

~  ~n'{n'  +  l)  ■ 

Finally,  from  equations  (42)  to  (44),  the  following  relations  are  valid: 


_  <(m,n) 

(45) 

n(m,n) 

“  -°n' 

(46) 
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6.  Plane  Interface 


The  boundary  conditions  at  the  plane  interface  are  the  most  difficult  to 
satisfy  because  all  the  field  expansions  are  in  terms  of  vector  spherical 
harmonics.  In  general  terms,  we  can  introduce  another  transition  matrix 
which,  for  our  purposes,  transforms  any  set  of  expanding  spher¬ 
ical  waves  centered  on  coordinate  system  {xi,yi,zi)  incident  on  the  plane 
interface: 

EElr  =  (fcin)  +  N®  (fein) ,  (47) 

to  its  scattering  response,  a  set  of  expanding  spherical  waves  centered  on 
coordinate  system  {x2,y2,  Z2)  that  are  reflected  by  the  plane  interface: 

cxD  n' 

■pplane  ^  L{1)  r/H  i  „(2)  rrl2  1  ^(3)  /,  n 

n'=l  ra'=—n' 

rnnrn' n'  ^mn^mnm'n'  ]  (^l'^2)  ,  (48) 

that  is,  any  field  incident  on  the  plane  interface  expressed  by  equation  (47) 
is  scattered  by  a  field  expressed  by  equation  (48).  For  a  plane  interface,  all 
elements  =  ^mm'U^^nm'n'  because  of  the  system  symmetry. 

We  first  treat  the  case  of  a  perfectly  conducting  plane  interface.  At  the  inter¬ 
face  (r2  =  ri,  6*2  =  vr— 01,  02  =  0i),  the  tangential  components  of  the  electric 
and  magnetic  fields  must  be  continuous.  We  can  achieve  this  by  making  the 
scattered  field  the  mirror  image  of  the  incident  field.  Using  equations  (1), 
(2),  (16),  and  (17),  we  can  easily  show  that  the  boundary  conditions  are  sat¬ 
isfied  if  the  scattered  field  is  of  the  form 

e£”  =  (/W  (-ir+"*+‘M(-ir+’"NW  (fcir,)  (49) 

or 

ULm'n'  =  (-1)"+™+^-  .  (50) 

This  is  equivalent  to  the  following  relationship  between  the  scattering  and 
interaction  coefficients 

=  (51) 

and  equation  (36)  can  now  be  written  as 
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Y  1  «mn  +  (“ 

n=l  m=—n  K  n'=l 

{oo 

<■£  +  E  (-!)"'*"• 

n'=l 

Equation  (52)  represents  the  final  solution  to  the  scattering  problem  of  a 
sphere  near  a  perfectly  conducting  interface.  The  scattering  coefficients  bmh 
are  written  in  terms  of  a  set  of  known  quantities:  translation  coefficients 
^  and  \  incident  field  coefficients  Umn,  and  Mie  coefficients  fl- 
By  truncating  the  infinite  series,  we  can  find  a  solution  through  matrix  in¬ 
version  techniques.  The  number  of  terms  necessary  for  convergence  de¬ 
pends  on  the  required  accuracy.  We  can  obtain  adequate  results  by  carrying 
the  summation  in  n  to 

N  =  Xsph  +  +  2  (53) 

terms,  where  Xgph  is  the  size  parameter  of  the  sphere.  Note  that  this  formula 
is  independent  of  the  separation  distance  d. 

Many  real-world  interfaces  are  not  perfectly  conducting.  To  solve  this  more 
general  problem,  we  must  either  make  some  approximation  or  resort  to  nu¬ 
merical  techniques.  The  simplest  approximation  one  can  make  is  to  assume 
that,  for  the  purposes  of  the  interaction  only,  the  image  fields  are  propor¬ 
tional  to  the  scattered  fields,  i.e., 

(0)  .  (54) 

The  factor  of  proportionality  chosen  is  the  Fresnel  coefficient  at  normal  in¬ 
cidence.  This  means  that  the  scattered  fields  from  the  sphere  that  reflect  off 
the  plane  interface  and  strike  the  sphere  again,  strike  the  plane  interface 
at  normal  incidence.  Geometrically,  there  is  some  justification  to  this  argu¬ 
ment:  rays  emanating  from  the  center  of  the  sphere  reflecting  off  the  plane 
interface  and  striking  the  sphere  can  have  a  maximum  incident  angle  on  the 
plane  interface  of  30° ,  and  the  Fresnel  reflection  coefficients  remain  nearly 
constant  at  near-normal  incident  angles.  Equation  (54)  is  equivalent  to  the 
following  relationship  between  the  scattering  and  interaction  coefficients 

=  -  (-1)“^’  Rji  (0)  6^  (55) 


,  Nn'+m  r,  (1)  A(m,n')  _  A2)  D(m,n')l  1  rl 

'^mn'^n  '^mn'^n  (  Jn 


l(1)  Tj{m,n')  _  7  (2)  .(m,n')l  1  r2 

^mn'^n  (  Jn  ■ 


(52) 


15 


and  equation  (36)  can  now  be  written  as 


oo  n 


n=l  m=—n 

j  o!.;i  +  E  (-1)”'+’”  (0)  .Af’”''  -  C.-Bi’"'"'’]  I  ft. 

I  n'=l  J 

+  ( <•£  +  E  (-1)”'*”  Bjf  (0)  I  fl 


n'=l 


An  exact  solution  for  the  scattering  from  a  sphere  located  near  an  arbitrary 
half  space  requires  considerably  more  effort;  however,  we  can  find  one  us¬ 
ing  the  derivation  above.  The  only  requirement  is  finding  the  matrix  ele¬ 
ments  U^^nm'n' '  ^7  expanding  the  individual  vector 

spherical  harmonics  given  by  equation  (47)  as  an  angular  spectrum  of  plane 
waves.  The  reflected  field  is  also  expressed  as  an  angular  spectrum  of  plane 
waves  and  the  angle-dependent  amplitudes  of  the  reflected  field  are  equal 
to  those  of  the  incident  field,  mulitplied  by  the  appropriate  Fresnel  reflec¬ 
tion  coefficient.  Finding  the  matrix  elements  requires  expressing 

the  reflected  angular  spectrum  in  terms  of  the  vector  spherical  harmonics. 
This  final  step  is  a  minor  variation  on  equations  (8)  and  (10).  Note  that  with 
little  modification,  the  equations  of  this  report  can  also  be  modified  to  ac¬ 
commodate  irregularly  shaped  particles  and  roughness  on  the  "plane"  in¬ 
terface.  The  only  knowledge  required  is  the  appropriate  T-matrix  of  these 
constituents. 
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7.  Scattered  Field 


For  many  applications,  only  the  far-field  region  is  of  interest,  where 

hW  (tr)  =  (57) 

ikr 

The  scattered  far  fields  can  be  expressed  in  terms  of  a  scattering  amplitude 
matrix  that  contains  all  the  polarization  information  of  the  scattered  light, 

(  \  ^  exp  [ifcin]  [52  *53  1  /  ^i™  \ 

I  EL  )  L  ^4  5,  J  I,  J  ’ 

where  ^tm  are  the  amplitudes  of  incident  TE  and  TM  plane 

waves,  and  EL  arid  EL  are  the  61  and  L  components  of  the  scattered  elec¬ 
tric  field.  The  total  scattered  field  is  composed  of  the  linear  superposition  of 
the  scattered  field  from  the  sphere  Egca  and  from  its  image  Eima-  We  have 
provided  only  an  explicit  expression  for  Eima  when  the  interface  is  per¬ 
fectly  conducting  (for  the  dielectric  half  space,  the  approximate  coefficients 
given  are  only  valid  for  the  purpose  of  the  interaction).  However,  the  con¬ 
tribution  of  the  image  field  in  the  far  field  is  equivalent  to  the  image  of  the 
scattered  field  multiplied  by  the  Fresnel  reflection  coefficient  i?TE  L  —  Oi  ) 
or  Rtm  (tt  —  0i),  and  a  phase  shift  exp  {—2ikdcos9i).  We  can  express  the 
image  of  the  scattered  field  in  terms  of  the  scattered  field  using  equations 
(16)  and  (17).  Separating  the  individual  vector  components  of  equation  (22) 
yields  the  following: 


=  (-i)”exp(im(?!>i) 


n=l  m=—n 


X  I  1  -|-  i?™  (tt  —  61)  (— 1)”+”^  exp  {—2ikdcos  9i)  (cos  9i) 

+  l-i?™(7r-6'i)(-l)”+™exp(-2i/c(icos6'i)6W,™fmn(cos6»i) 


=  £  £  (-i)”+^exp(im(^i) 


n=l  m=—n 


X  I  1  -I-  i?™  (tt  —  9i)  {—l)'^~^'^exp{—2ikdcos9i)  (cos0i) 


+  [l-i?™(7r-0i)(-l) 


exp  {—2ikdcos  9i)  Lin. 


nm  '  Tnn 


(cos  01 )  I , 
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(61) 


oo  n 

*^3  =  XI  XI  (-i)”’^^exp(im(/>i) 

n=l  m=—n 

x|  1  +  R™  {tt  —  6i)  {—1)'^'^'^  exp{—2ikdcos9i)  {cos  6i) 

+  1 -i?™(7r-6'i)(-l)"'’^™exp(-2i/c(icos6'i)6^^j,™fmn(cos6'i) 

OO  n 

54  =  X  X  (-i)”exp(im(/>i)  (62) 

n=l  m=—n 

X  I  1  +  i?™  (vr  —  9i)  (—1)"'+™'  exp  (— 2iA:dcos  9i)  (cos  9i) 

+  1  - -R™  (tt  -  0i)  (-l)"^’''”^exp  (-2iA:(icos6'i)  (cos6*i) 

The  superscripts  on  the  scattering  coefficients  bn)^^  and  bn)J^  are  the  scat¬ 
tering  field  coefficients  when  the  system  is  illuminated  by  TE  and  TM  plane 
waves.  These  expressions  represent  the  far-field  solution  to  the  problem  of 
a  sphere  near  a  plane  interface.  Scattering  amplitudes  are  convenient,  the¬ 
oretically,  because  they  can  be  manipulated  to  produce  measureable  quan¬ 
tities  like  scattering  intensities  and  polarization  states.  Figure  2  shows  the 
scatter  in  the  plane  of  incidence  {(pi  =  0°)  of  r  =  0.55  pm  spheres  rest¬ 
ing  on  a  substrate  illuminated  at  9o  =  30°  and  A  =  0.6328  pm.  This  fig¬ 
ure  demonstrates  the  sensitivity  of  the  scatter  to  system  parameters.  The 
TM  polarization  intensities  are  especially  sensitive  to  a  change  in  the  half¬ 
space  complex  refractive  index  (m2),  since  the  TM  Fresnel  reflection  am¬ 
plitudes  are  more  sensitive  than  the  TE  amplitudes.  Many  sizing  routines 
are  based  on  the  positions  and  spacing  of  the  maxima  and  minima  intensi¬ 
ties.  This  graph  demonstrates  that  care  must  be  taken  when  applying  such 
techniques  to  particles  located  on  substrates.  The  spacings  remain  fairly 
constant  for  the  TE  illumination  of  figure  2,  but  change  drastically  for  TM 
illumination. 
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Figure  2.  Light  scattering  intensities  calculated  for  r  —  0.55  spheres  on  a  substrate 
illuminated  at  A  =  0.6328  /im,  6o  =  30°,  4>i  =0°,  and  d  —  0.55  /rm.  Results  are  shown  when 
system  is  composed  entirely  of  gold  (rriaph  =  0.226  +  3.32i)  and  when  sysfem  is  perfectly 
conducting. 
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8.  Conclusion 


We  derived  a  solution  for  the  scattering  from  a  spherical  particle  near  a 
perfectly  conducting  plane  interface.  The  analytical  solution  is  in  terms  of 
a  vector  spherical  harmonic  expansion.  The  equations  can  be  modified  to 
accommodate  an  irregular  particle  near  an  irregular  interface  if  the  scat¬ 
tering  T-matrix  is  known  for  these  isolated  subsystems.  Additionally,  the 
equations  may  also  accommodate  a  nonperfectly  conducting  half  space.  We 
provided  an  approximate  solution  for  this  latter  case,  based  on  a  normal- 
incident  Fresnel  reflection  of  the  interaction  field. 
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